First-principles study of the electronic structure and magnetism of CalrOa 
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I study the electronic structure and magnetism of postperovskite CalrOg using first-principles 
calculations. The density functional calculations within the local density approximation without 
the combined effect of spin-orbit coupling and on-site Coulomb repulsion show the system to be 
metallic, which is in disagreement with the recent experimental evidences that show CalrOa to be 
an antiferromagnetic Mott insulator in the Jeff =1/2 state. However, when spin-orbit coupling is 
taken into account, the Ir t2g bands split into fully filled J^a = 3/2 bands and half-filled J^s = 1/2 
bands. I find that spin-orbit coupling along with a modest on-site Coulomb repulsion opens a gap 
leading to a Mott insulating state. The ordering is antiferromagnetic along the c axis with total 
moments aligned antiparallel along the c axis and canted along the b axis. 

PACS numbers: 71.30.+h, 75.25.Dk, 75.50.Ee 



I. INTRODUCTION 



Transition-metal oxides (TMOs) in perovskitc and 
related structures exhibit myriad interesting proper- 
ties. These include unconventional superconductivity 
in cuprateSfi colossal magnetoresistance in manganites?^ 
and fcrroelectricity in Pb(Zr,Ti)034 In these materials, 
the transition-metal ion is situated inside an oxygen octa- 
hedral cage, which may be arranged in a corner- or edge- 
shared manner. The interesting properties of TMOs arise 
because of the competition between the crystal-field split- 
ting (which arises because of covalency between transi- 
tion metal d and oxygen p states), on-site Coulomb repul- 
sion [/, Hund's coupling, spin-orbit (SO) coupling due to 
orbital degeneracy (leading to unquenched angular mo- 
ment in the ground state), and different spin-exchange 
pathways. 

In this regard, the discovery of a spin-orbital Mott 
state in Sr2lr04 by Kim et alA>^ is significant because 
it enables us to study the case where spin-orbit cou- 
pling and its interplay with the Coulomb repulsion is an 
important ingredient in determining the electronic and 
magnetic properties of the system. Sr2lr04 exists in a 
layered perovskitc structure. The Ir^"*" (5c?^) ions are sit- 
uated inside corner-shared O octahedral cages, which arc 
themselves arranged in a square lattice in the ah plane. 
As there arc an odd number of electrons per formula 
unit, one might expect this material to be a metal in 
the band picture. However, Sr2lr04 is experimentally a 
canted antiferromagnetic insulator As explained lucidly 
by Kim et al., a state with an effective total angular 
momentum Jeff —1/2 that has a complex wave func- 
tion is realized in Sr2lr04, which arises due to the com- 
bined influence of strong spin-orbit coupling and mod- 
erate on-site Coulomb repulsion A similar conclusion 
on Sr2lr04 has been reached through the study of a 
three-orbital Hubbard model with spin-orbit coupling^ 
and combined density functional theory and dynamical 
mean field theory (LDA+DMFT) calculations,-^ although 
Arita et al. suggest that Sr2lr04 is a Slater insulator 
based on their LDA-fDMFT study J. Interestingly, it was 



also found that spin-orbit coupling plays an important 
role in the electronic properties even for a AdP system 
such as Sr2Rh04 3iia 

It has been suggested that systems in the Jeff = 1/2 
state, depending on bond geometry, lead to interest- 
ing varieties of low-energy Hamiltonians, including the 
isotropic Heisenberg model and the highly anisotropic 
quantum compass or Kitaev models relevant for quantum 
computingjii Therefore, it is important to investigate 
materials that exhibit the Jeff —1/2 state in different 
structures in order to study the effect of different local en- 
vironments and spin-exchange pathways. Recently, Oh- 
gushi et al. reported resonant x-ray diffraction study of 
CalrOs that indicates this material also exhibits a Mott 
insulating Jeff =1/2 stateJ^ CalrOa exists in the post- 
perovskite structure with space group Cmcm as shown 
in Fig. [TJ The Ir^"*" (5d^) ions are situated inside the O 
octahedra, but these octahcdra share an edge along the c 
axis, unlike the case of Sr2lr04. Thus, CalrOa is another 
ideal material to investigate the interplay between spin- 
orbit coupling and on-site Coulomb repulsion that may 
help in understanding the unique properties that might 
be exhibited by Jeff =1/2 systems. 

Tsuchiya et al. have reported first-principles density 
functional calculations that show this material is a metal 
within the local density approximation^^ This is con- 
trary to the experimental evidence that shows this ma- 
terial is a Mott insulator that undergoes an antiferro- 
magnetic transition at Tn = 115 KJ^ A recent reso- 
nant x-ray diffraction study shows that the ordering is 
of stripe-type antiferromagnetism along the c axis, with 
total moments aligning parallel along the a axis and 
antiparallel along the c axisf^ The inverse susceptibil- 
ity 1/x deviates from the linear behavior at a temper- 
ature ~350 K that is considerably higher than Tn and 
a Curie- Weiss fit to x above 400 K gives a Curie- Weiss 
temperature of 3900 This indicates that the anti- 
ferromagnetic correlations arise much before the anti- 
ferromagnetic transition, and magnetic ordering is sup- 
pressed by low dimensionality or competing ordering in- 
teractions. Jang et al. have studied the electronic struc- 
ture of meta-stable perovskitc Cai-^^SrjjIrOs {x = 0, 0.5, 
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FIG. 1: (Color online) Crystal structure of CalrOs. The large 
(cyan) balls are Ca, small (red) balls are O, and the Ir atoms 
reside inside the (brown) octahedra. 




and 1) thin films using transport measurements, optical 
spectroscopy, and pseudopotential-based first-principles 
calculations They find that perovskite CalrOa thin 
films are semimetallic and near the metal-insulator tran- 
sition. Their calculations with spin-orbit coupling and 
on-site Coulomb repulsion U found that the spin-orbit 
couphng splits the Ir t2g states into Jeff = 3/2 and 1/2 
states, and U = 2.0 eV further splits the Jeff =1/2 states, 
although the valence and conduction bands still touch the 
Fermi level, resulting in a semimetallic state. 

The experimental evidences that have so far been ac- 
cumulated suggest that calculations that include the ef- 
fect of spin-orbit coupling and on-site Coulomb repulsion 
would be helpful in clarifying the electronic and magnetic 
properties of CalrOa. In this paper, I report the results 
of density functional calculations that show CalrOa is 
in a Mott insulating state that is induced by the com- 
bined effect of spin-orbit coupling and on-site Coulomb 
repulsion. This state arises out of spin-orbit split Ir t2g 
bands that get separated into lower lying fully filled and 
higher lying half-filled bands that have effective total an- 
gular momenta JeS = 3/2 and 1/2, respectively, in the 
strong spin-orbit coupling limit. The half-filled Jeff = 
1/2 bands are narrow, so even a modest on-site Coulomb 
repulsion induces a Mott insulating state that is topolog- 
ically different from the metallic state given by the lo- 
cal density approximation, without taking into account 
the spin-orbit coupling and on-site Coulomb repulsion. 
This is a Mott insulating state in the sense that a single- 
particle theory such as the density functional theory im- 
plemented using Kohn-Sham formalism cannot explain 
the insulating state, and an explicit treatment of on-site 
Coulomb repulsion is needed. The Mott insulating state 
thus obtained is antiferromagnetically ordered along the 
c axis with total moments aligned antiparallel along the 
c axis and canted along the b axis. 



II. APPROACH 

The purpose of this paper is to elucidate the role of 
spin-orbit coupling and on-site Coulomb repulsion on 



FIG. 2: (Color online) Non-spin-polarized LDA DOS of 
CalrOa (in states/eV). The projections are onto the respective 
muffin-tin spheres and are only indicative of the contribution 
to the total DOS. The Fermi energy is at eV. 



the electronic and magnetic properties of CalrOs us- 
ing density functional calculations. The calculations 
were performed within the local density approximation 
(LDA) using the general full-potential linearized aug- 
mented plane-wave method as implemented in the ELK 
software packageJ^ Muffin-tin radii of 2.2, 2.0, and 1.6 
a.u. for Ca, Ir, and O, respectively, were used. A 8 x 8 x 8 
/c-point grid was used to perform the Brillouin zone inte- 
gration, and the convergence of moments was checked on 
a 10 X 10 X 10 grid. The effect of spin-orbit coupling was 
treated using a second-variational scheme, and the fully 
localized limitil is used to take into account the double 
counting in LDA+C/ calculations. A value for the on-site 
Coulomb repulsion U = 2.75 cV (which gives a band gap 
close to the experimental value) was used unless other- 
wise mentioned. 

I used the experimental lattice parameters a = 3.145 
A, 5 = 9.855 A, and c = 7.293 A^ but relaxed the atomic 
positions. The calculated atomic positions Ca (0, 0.2498, 
0.25), Ir (0, 0, 0), 0(1) (0.5, 0.4253, 0.25), and 0(2) (0.5, 
0.1230, 0.0485) agree well with the experimental values 
Ca (0, 0.2498, 0.25), Ir (0, 0, 0), 0(1) (0.5, 0.4331, 0.25), 
and 0(2) (0.5, 0.1296, 0.0553). The results presented in 
this paper are for the relaxed atomic positions, but I also 
performed calculations with experimental atomic posi- 
tions and came to the same physical conclusions. There 
are two formula units per primitive unit cell in the Cmcm 
structure. The Ir^"*" ions make a two-dimensional rectan- 
gular lattice in the ac plane (not shown) and the Ir-0 
layer is separated by a layer of Ca atoms along the b 
axis. For the relaxed atomic positions, the O octahe- 
dra are tilted by an angle of 22°. The O octahedra are 
slightly compressed along the corner-shared O direction 
(from left to right in the left figure of Fig. [1]) with a 
bond-length ratio of 0.97. The Ir-0 distances along the 
corner-shared c axis and edge-shared a axis are 1.97 and 
2.02 A, respectively. 
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III. RESULTS 

Let us first consider the non-spin-polarized (NSP) LDA 
calculations. Even though they are inadequate to de- 
scribe the ground-state properties, these calculations give 
a decent description of the band structure of CalrOa that 
provide a playground for the interplay between SO cou- 
pling and U. The electron density of states (DOS) and 
band structure for this case are shown in Figs. [2] and 
[3l^a), respectively. Most of the bands along the b direc- 
tion (T-Y) have low dispersion, which suggests that the 
physics related to two dimensionality might be relevant 
in this system. There are 24 bands between —9.0 and 1.0 
eV. The 18 bands between —9.0 and —2.3 eV (not shown) 
have a dominant O p character and thus derive from the 
p orbitals of the six O atoms in the unit cell. These bands 
also show Ir d character, which implies significant cova- 
lency between the O p and Ir d states as the unoccupied 
Ir d states above Fermi level also contain some admixture 
of O p states. There is a very small gap of ~0.05 eV at 
—2.36 eV; beyond which lie six bands with a dominant Ir 
d character. These are the Ir t2g states that are formally 
antibonding, and these bands correspondingly show some 
O p contribution. A gap of '^1.5 eV separates these t2g 
states from a group of four bands that have a mostly Ir 
d character, which are the Ir Cg states. The Ir Cg also 
have some O p character due to Ir d-O p covalency. The 
Ir 5d states are quite delocalized and the edge-sharing 
compressed IrOe octahedra are rotated by 22°, and this 
leads to some hybridization the between Ir t2g and Cg 
levels. The Ca and Ir s states are high above the Fermi 
level, and within an ionic limit the electronic structure is 
consistent with the ionic states Ir''"'" and 0^~, although 
there is significant deviation from this because of Ir d-0 
p hybridization. The two Ir^"*" ions nominally have five 
electrons each in their d orbitals. As a result, the six Ir 
t2g bands are not fully filled, and the system is a metal 
within LDA with a t2g hole on each Ir^"^ ion. 

The metallic state obtained within LDA due to incom- 
plete filling of Ir t2g states is contrary to the experimental 
evidence that indicates CalrOs is a Mott insulator. This 
suggests that spin-orbit coupling and/or on-site Coulomb 
repulsion play crucial roles in the electronic and magnetic 
properties of CalrOa. Let us now consider the effect of 
SO coupling and U on the electronic structure of CalrOa. 
The NSP LDA, LDA+C/, LDA-hSO, and LDA-hS0-hJ7 
(with U = 2.75 eV) Ir t2g bands are shown in Fig. [3l 
Let us first note that a value for U of 2.75 eV without 
the spin-orbit coupling has very little effect on the band 
structure (I did the calculation with [/ up to 5 eV with- 
out getting an insulating state). This is not surprising as 
the Ir <2g manifold is spread over a bandwidth of ^--^2.8 
eV, and it would require a substantially larger U to open 
up a gap. However, turning on spin-orbit coupling makes 
a significant difference in the electronic structure. The 
system is non-magnetic, so the bands are spin degener- 
ate as they are not exchange split. However, the spin- 
oribt coupling splits the manifold of six spin-degenerate 
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FIG. 3: NSP LDA, LDA+U, LDA+SO and LDA-f-SO+f/ 
band structures, respectively from top to bottom, of CalrOs. 
The band structures are plotted along the path F (0,0,0) — !■ Y 
(0,i,0) ^ T (0, i, i) ^ r (0,0,0) ^ S (i, i 0) ^ R (i, i, 
i) ^ r (0,0,0) Z (0,0,i). Here, U = 2.75 eV is used. The 
bands are exchange split only for the case of LDA+SO-l-17. 



Ir t2g bands into a lower-lying group of four and a higher- 
lying group of two spin-degenerate bands. In the limit 
of strong spin-orbit coupling, the lower and higher sets 
of bands within the t2g manifold would correspond to 
effective total angular momenta Jcs of 3/2 and 1/2, re- 
spectively. This is similar to the case of Sr2lr04 where 
the spin-orbit coupling splits the Ir t2g bands into a lower 
lying quartet of J^s = 3/2 and a higher lying doublet of 
J^ff =1/2 bandsi^ In the case of CalrOa, the J^s =1/2 
bands are narrow with a width of ^1 eV. These are sep- 
arated from the Jctf = 3/2 bands by ~0.15 eV, although 
the gap is indirect. The ten Ir d electrons in the unit cell 
completely fill the Jos = 3/2 bands, while the Jcs =1/2 
bands are only half filled. As a result, the system is still 
metallic. 

Even though on-site Coulomb repulsion U and spin- 
orbit coupling acting alone do not make the system an 
insulator, it is likely that their combined effect can in- 
duce a Mott insulating state by splitting the narrow Jeff 
= 1/2 bands. The LDA+SO+C/ calculations with U = 
2.75 eV reveal that this scenario is realized in CalrOa. 
As shown in Fig. [3]^d), the U in the presence of SO cou- 
pling makes only minor modifications to the Jeff = 3/2 
bands. The Jeff = 3/2 bands get exchange split and a 
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TABLE I: The {L) and {S) expectation values computed over Ir mufRn-tin spheres and the band gap -Bgap (eV) for some values 
of on-site Coulomb repulsion U (eV) and Hund's coupling J (eV). The moments are in units of Bohr magneton. 



Site 


{L) 


{S) 




Ir(l) 
Ir(l') 


(0.00,0.06, -0.28) 
(0.00,0.06,0.28) 


(0.00,0.03,-0.19) 
(0.00,0.03,0.19) 


U = 2.75, J = 0.0 

^'gap =0.33 


Ir(l) 
Ir(l') 


(0.00,0.06, -0.27) 
(0.00,0.06,0.27) 


(0.00,0.03,-0.18) 
(0.00,0.03,0.18) 


U = 2.75, J = 0.3 
£gap =0.30 


Ir(l) 
Ir(l') 


(0.00,0.04,-0.18) 
(0.00,0.04, 0.18) 


(0.00,0.02,-0.13) 
(0.00,0.02,0.13) 


(7 = 2, J = 0.0 
£gap = 0.03 



degeneracy at the point T (0, 0.5, 0.5) is lifted, but oth- 
erwise the bandwidth and topology of the bands do not 
change substantially. However, the half- filled Jeff =1/2 
bands, in addition to being exchange split by ~0.1 eV, 
are split into the upper (UHB) and lower (LHB) Hub- 
bard bands, yielding a Mott insulating state. The fully 
occupied Jeff = 1/2 LHB has a small bandwidth of ^0.5 
eV and is separated from the unoccupied UHB by a gap 
of -0.33 eV for U = 2.75 eV (the gap is -0.03 eV for 
t/ = 2 eV). This agrees well with a band gap of 0.34 
eV obtained experimentally] "'^^'^^ I also performed calcu- 
lations with U = 1.0, 1.5, 2.0, and 2.5 eV. The system 
is nictalhc within LDA-I-SO-I-C/ for U up to 1.5 eV, but 
it becomes an insulator by J7 = 2.0 eV. To see the effect 
of the Hund couphng J, I did LT)A+SO+U calculations 
with U = 2.75 eV and J = 0.1, 0.2, and 0.3 eV. I find 
that these values of Hund coupling J do not change the 
qualitative picture — the Ir t2g levels are still split into Jeff 
= 3/2 and 1/2 bands and the Jeff =1/2 bands are fur- 
ther split into fully occupied LHB and unoccupied UHB. 
The inclusion of Hund coupling mainly reduces the band 
gap (for J = 0.3 cV, the band gap is 0.30 eV) and the 
magnetic moment. 

The LDA+SO+C/ calculations give an antiferromag- 
netic ground state for CalrOs along the c axis with total 
moments aligning antiparallcl along the c axis. The or- 
bital and spin moments are parallel to each other along 
the c axis and ferromagnetically canted along the b axis. 
The angular and spin expectation values computed over 
the two Ir muffin-tin spheres Ir(l) (0.0, 0.0, 0.0) and Ir(l') 
(0.0, 0.0, 0.5) for different U and J values are given in 
TableHl For U = 2.75 eV, the total moment is 0.67 ^s/Ir 
with an orbital moment of 0.29 fj,B and a spin moment 
of 0.38 (= 2|(iS')|). The canting angle is 10°, approxi- 
mately half the octahedral tilting angle of 22° and twice 
the value of 4° reported in Ref. [T^- The calculated val- 
ues differ considerably from what is expected for the ideal 
Jeff =1/2 state. In the ionic limit, one expects an orbital 
moment of 0.67 /is and a spin moment of 0.33 fiB foi' a- 
Jeff =1/2 statci^ In contrast, I obtain an orbital moment 
that is lower than the spin moment. The reason for this 
deviation from the Jeff = 1/2 may be the compression 
and tilting of the IrOg octahedra, in addition to the co- 
valency between Ir d and O p states. The compression of 
the IrOg octahedra will quench the orbital moment as the 
degeneracy between t2g states arc lifted. Also, the tilting 



of the IrOg octahedra causes the Cg bands to get mixed 
with the t2g states. These two effects should reduce the 
orbital moment but might enhance the spin contribution. 
It is interesting to note that Sr2lr04 also has distortion 
of the IrOg octahedra with a bond-length ratio of 1.04 
and a tilting angle of 11°^ and it has a calculated or- 
bital moment of 0.26 ^b and spin moment of 0.10 ^b^^ 
As the tilting angle in CalrOa is twice that of Sr2lr04, 
it might be reasonable to expect that CalrOs deviates 
further from the ideal Jeff =1/2 state due to the mixing 
of the Cg states. 



IV. CONCLUSIONS 

In summary, the electronic structure and magnetic 
properties of CalrOs has been studied using first- 
principles calculations. The system is metallic within 
the LDA because the Ir t2g states are incompletely filled. 
Modest values of on-site Coulomb repulsion alone have 
very little effect on the LDA electronic structure as the 
Ir t2g states have a broad bandwidth. The introduction 
of spin-orbit coupling splits the Ir t2g states into fully 
filled Jeff = 3/2 bands and half-filled Jeff =1/2 bands. 
The half- filled bands have a small bandwidth of —1.0 eV, 
which is split into fully filled lower and unfilled upper 
Hubbard bands by a modest value of on-site Coulomb 
repulsion. This topologically non-trivial modification of 
the LDA electronic structure due to the combined effects 
of spin-orbit coupling and on-site Coulomb repulsion re- 
sults in a Mott insulating state that is antiferromagnet- 
ically ordered along c axis with total moments aligning 
antiparallcl along the c axis and canted along the b axis. 
For U = 2.75 cV, the total magnetic moment is 0.67 
fiB with an orbital contribution of 0.28 /is and a spin 
contribution of 0.38 fiB- These values differ from what is 
expected for the ideal Jeff =1/2 state, and this deviation 
might be explained by the mixing of Jeff =1/2 bands with 
Ir eg bands due to the tilting of IrOg octahedra. There 
has been great interest in finding different materials with 
unique magnetic properties. The results presented here 
give strong support to the claim made by Ohgushi et al. 
in Ref. [T^ that CalrOa has a unique spin-orbit integrated 
magnetic ground state. 



5 



V. ACKNOWLEDGEMENTS L. Boeri for encouragement, helpful discussions, and sug- 

gestions in improving the manuscript, and to A. Avella 
I am grateful to G. Jackeh, D. I. Khomskii, D. J. Singh, for critical reading of the manuscript, 
and I. I. Mazin for helpful discussions. I am thankful to 



^ J. G. Bednorz and K. A. MuUer, Z. Phys. B 64, 189 (1986). 
^ R. von Helmolt, J. Wecker, B. Holzapfel, L. Schultz, and 

K. Samwer, Phys. Rev. Lett. 71, 2331 (1993). 
^ B. Jaffe, W. J. Cook, and J. Jaffe, Piezoelectric Ceramics 

(Academic, London, 1971). 
" B. J. Kim, H. Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, C. 

S. Leem, J. Yu, T. W. Noh, C. Kim, S.-J. Oh, J.-H. Park, 

V. Durairaj, G. Cao, and E. Rotenberg, Phys. Rev. Lett 

101, 076402 (2008). 
^ B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, 

and H. Takagi, Science 323, 1329 (2009). 
® H. Watanabe, T. Shirakawa, S. Yunoki, Phys. Rev. Lett. 

105, 216410 (2010). 

C. Martins, M. Aichhorn, L. Vaugier, and S. Biermann, 

Phys. Rev. Lett. 107, 266404 (2011). 
^ R. Arita, J. Kunes, A. V. Kozhevnikov, A. G. Eguiluz, and 

M. Imada, larXiv: 1107.0835] 
^ M. W. Haverkort, L S. Elfimov, L. H. Tjeng, G. A. 

Sawatzky, and A. Damascelh, Phys. Rev. Lett. 101, 026406 

(2008) . ' 

" G.-Q. Liu, V. N. Antonov, O. Jepsen, and O. K. Andersen, 

Phys. Rev. Lett. 101, 026408 (2008). 
^ G. Jackeh and G. KhaliuUin, Phys. Rev. Lett. 102, 017205 

(2009) . 



K. Ohgushi, J. -I. Yamaura, H. Ohsumi, K. Sugimoto, S. 
Takesh ita, A. Tokuda, H . Takagi, M. Takata, and T.-H. 
Arima, larXiv:1108.4523] 
^ T. Tsuchiya and J. Tsuchiya, Pliys. Rev. B 76, 144119 
(2007). 

" K. Ohgushi, H. Gotou, T. Yagi, Y. Kiuchi, F. Sakai, and 

Y. Ueda, Phys. Rev. B 74, 241104(R) (2006). 

S. Y. Jang, H. Kim, S. J. Moon, W. S. Choi, B. C. Jeon, J. 

Yu, and T. W. Noh, J. Phys.: Condens. Matter 22, 485602 

(2010). 

^ http: / /elk. sourcef orge.net] 

V. L Anisimov, I. V. Solovyev, M. A. Korotin, M. T. 

Czyzyk, and G. A. Sawatzky, Phys. Rev. B 48, 16929 

(1993). 

* R. von Fritz and B. Dietrich, Z. Anorg. AUg. Chem. 336, 
17 (1965V 

^ In Ref. Ujj . the resistivity versus temperature curve is fit to 
the formula p{T) = A exp(A/r) and a value of A = 0.17 
eV is obtained. The band gap should be twice the value of 
A, i.e. 0.34 eV. 

'° M. K. Crawford, M. A. Subramanian, R. L. Harlow, J. A. 
Fernandez-Baca, Z. R. Wang, and D. C. Johnston, Phys. 
Rev. B 49, 9198 (1994). 



